Autoimmunity Against Surfactant Protein B Is Associated with Pneumonitis During Checkpoint Blockade.

RATIONALE: Immune checkpoint inhibitor-related pneumonitis is a serious autoimmune event affecting up to 20% of patients with non-small cell lung cancer, yet the factors underpinning its development in some patients and not others are poorly understood. OBJECTIVES: To investigate the role of autoantibodies and autoreactive T cells against surfactant-related proteins in the development of pneumonitis. METHODS: The study cohort consisted of non-small cell lung cancer patients who gave blood samples before and during immune checkpoint inhibitor treatment. Serum was used for proteomics analyses and to detect autoantibodies present during pneumonitis. T cell stimulation assays and single-cell RNA sequencing were performed to investigate the specificity and functionality of peripheral autoreactive T cells. The findings were confirmed in a validation cohort comprising patients with non-small cell lung cancer and patients with melanoma. MEASUREMENTS AND MAIN RESULTS: Across both cohorts, patients who developed pneumonitis had higher pre-treatment levels of immunoglobulin G autoantibodies targeting surfactant protein-B. At the onset of pneumonitis, these patients also exhibited higher frequencies of CD4+ interferon-gamma-positive surfactant protein B-specific T cells, and expanding T cell clonotypes recognizing this protein, accompanied by a pro-inflammatory serum proteomic profile. CONCLUSIONS: Our data suggest that the co-occurrence of surfactant protein-B-specific immunoglobulin G autoantibodies and CD4+ T cells is associated with the development of pneumonitis during ICI therapy. Pre-treatment levels of these antibodies may represent a potential biomarker for elevated risk of developing pneumonitis and on-treatment levels may provide a diagnostic aid.

Disruption of a Plasmodium falciparum patatin-like phospholipase delays male gametocyte exflagellation.

An essential process in transmission of the malaria parasite to the Anopheles vector is the conversion of mature gametocytes into gametes within the mosquito gut, where they egress from the red blood cell (RBC). During egress, male gametocytes undergo exflagellation, leading to the formation of eight haploid motile microgametes, while female gametes retain their spherical shape. Gametocyte egress depends on sequential disruption of the parasitophorous vacuole membrane and the host cell membrane. In other life cycle stages of the malaria parasite, phospholipases have been implicated in membrane disruption processes during egress, however their importance for gametocyte egress is relatively unknown. Here, we performed comprehensive functional analyses of six putative phospholipases for their role during development and egress of Plasmodium falciparum gametocytes. We localize two of them, the prodrug activation and resistance esterase (PF3D7_0709700) and the lysophospholipase 1 (PF3D7_1476700), to the parasite plasma membrane. Subsequently, we show that disruption of most of the studied phospholipase genes does neither affect gametocyte development nor egress. The exception is the putative patatin-like phospholipase 3 (PF3D7_0924000), whose gene deletion leads to a delay in male gametocyte exflagellation, indicating an important, albeit not essential, role of this enzyme in male gametogenesis.

The history and organization of the Workshop on Population and Speciation Genomics.

With the advent of high-throughput genome sequencing, bioinformatics training has become essential for research in evolutionary biology and related fields. However, individual research groups are often not in the position to teach students about the most up-to-date methodology in the field. To fill this gap, extended bioinformatics courses have been developed by various institutions and provide intense training over the course of two or more weeks. Here, we describe our experience with the organization of a course in one of the longest-running extended bioinformatics series of workshops, the Evomics Workshop on Population and Speciation Genomics that takes place biennially in the UNESCO world heritage town of Ceský Krumlov, Czech Republic. We list the key ingredients that make this workshop successful in our view, explain the routine for workshop organization that we have optimized over the years, and describe the most important lessons that we have learned from it. We report the results of a survey conducted among past workshop participants that quantifies measures of effective teaching and provide examples of how the workshop setting has led to the cross-fertilisation of ideas and ultimately scientific progress. We expect that our account may be useful for other groups aiming to set up their own extended bioinformatics courses.

The time is ripe for functional genomics: Can epigenetic changes mediate reproductive timing?

Populations are under strong selection to match reproductive timing with favourable environmental conditions. This becomes particularly important and challenging with increasing interannual environmental variability. Adjusting reproductive timing requires the ability to sense and interpret relevant environmental cues, while responding flexibly to their interannual variation. For instance, in seasonal species, reproductive timing is often dependent on photoperiod and temperature. Although many genes influencing the timing of reproduction have been identified, far less attention has been paid to the gene-regulatory cascades orchestrating these complex gene-environment interactions. In a From the Cover article in this issue of Molecular Ecology, Lindner, Laine, et al. (2021) addressed this knowledge gap by investigating the role of DNA methylation in mediating reproductive timing in the seasonally breeding great tit (Parus major). Using a clever blood sampling design, they investigated genome-wide DNA methylation changes following individual female birds across multiple reproductive stages. This approach revealed 10 candidate genes with a strong correlation between promoter methylation and reproductive status. Some of these genes are known to be involved in reproductive timing (e.g., MYLK-like or NR5A1), yet for others this function was previously unknown (Figure 1). Interestingly, NR5A1 is a key transcription factor, which may affect other genes that are part of the same regulatory network. The findings of Lindner, Laine, et al. (2021) provide a strong case for studying DNA methylation to uncover how gene-environment interactions influence important life-history traits, such as reproductive timing.

Individual variability and versatility in an eco-evolutionary model of avian migration.

Seasonal migration is a complex and variable behaviour with the potential to promote reproductive isolation. In Eurasian blackcaps (Sylvia atricapilla), a migratory divide in central Europe separating populations with southwest (SW) and southeast (SE) autumn routes may facilitate isolation, and individuals using new wintering areas in Britain show divergence from Mediterranean winterers. We tracked 100 blackcaps in the wild to characterize these strategies. Blackcaps to the west and east of the divide used predominantly SW and SE directions, respectively, but close to the contact zone many individuals took intermediate (S) routes. At 14.0° E, we documented a sharp transition from SW to SE migratory directions across only 27 (10-86) km, implying a strong selection gradient across the divide. Blackcaps wintering in Britain took northwesterly migration routes from continental European breeding grounds. They originated from a surprisingly extensive area, spanning 2000 km of the breeding range. British winterers bred in sympatry with SW-bound migrants but arrived 9.8 days earlier on the breeding grounds, suggesting some potential for assortative mating by timing. Overall, our data reveal complex variation in songbird migration and suggest that selection can maintain variation in migration direction across short distances while enabling the spread of a novel strategy across a wide range.

Experimental Parasite Infection Causes Genome-Wide Changes in DNA Methylation.

Parasites are arguably among the strongest drivers of natural selection, constraining hosts to evolve resistance and tolerance mechanisms. Although, the genetic basis of adaptation to parasite infection has been widely studied, little is known about how epigenetic changes contribute to parasite resistance and eventually, adaptation. Here, we investigated the role of host DNA methylation modifications to respond to parasite infections. In a controlled infection experiment, we used the three-spined stickleback fish, a model species for host-parasite studies, and their nematode parasite Camallanus lacustris. We showed that the levels of DNA methylation are higher in infected fish. Results furthermore suggest correlations between DNA methylation and shifts in key fitness and immune traits between infected and control fish, including respiratory burst and functional trans-generational traits such as the concentration of motile sperm. We revealed that genes associated with metabolic, developmental, and regulatory processes (cell death and apoptosis) were differentially methylated between infected and control fish. Interestingly, genes such as the neuropeptide FF receptor 2 and the integrin alpha 1 as well as molecular pathways including the Th1 and Th2 cell differentiation were hypermethylated in infected fish, suggesting parasite-mediated repression mechanisms of immune responses. Altogether, we demonstrate that parasite infection contributes to genome-wide DNA methylation modifications. Our study brings novel insights into the evolution of vertebrate immunity and suggests that epigenetic mechanisms are complementary to genetic responses against parasite-mediated selection.

Two different epigenetic information channels in wild three-spined sticklebacks are involved in salinity adaptation.

Epigenetic inheritance has been proposed to contribute to adaptation and acclimation via two information channels: (i) inducible epigenetic marks that enable transgenerational plasticity and (ii) noninducible epigenetic marks resulting from random epimutations shaped by selection. We studied both postulated channels by sequencing methylomes and genomes of Baltic three-spined sticklebacks (Gasterosteus aculeatus) along a salinity cline. Wild populations differing in salinity tolerance revealed differential methylation (pop-DMS) at genes enriched for osmoregulatory processes. A two-generation experiment demonstrated that 62% of these pop-DMS were noninducible by salinity manipulation, suggesting that they are the result of either direct selection or associated genomic divergence at cis- or trans-regulatory sites. Two-thirds of the remaining inducible pop-DMS increased in similarity to patterns detected in wild populations from corresponding salinities. The level of similarity accentuated over consecutive generations, indicating a mechanism of transgenerational plasticity. While we can attribute natural DNA methylation patterns to the two information channels, their interplay with genomic variation in salinity adaptation is still unresolved.

An exploration of the links between parasites, trophic ecology, morphology, and immunogenetics in the Lake Tanganyika cichlid radiation.

Differences in habitat and diet between species are often associated with morphological differences. Habitat and trophic adaptation have therefore been proposed as important drivers of speciation and adaptive radiation. Importantly, habitat and diet shifts likely impose changes in exposure to different parasites and infection risk. As strong selective agents influencing survival and mate choice, parasites might play an important role in host diversification. We explore this possibility for the adaptive radiation of Lake Tanganyika (LT) cichlids. We first compare metazoan macroparasites infection levels between cichlid tribes. We then describe the cichlids' genetic diversity at the major histocompatibility complex (MHC), which plays a key role in vertebrate immunity. Finally, we evaluate to what extent trophic ecology and morphology explain variation in infection levels and MHC, accounting for phylogenetic relationships. We show that different cichlid tribes in LT feature partially non-overlapping parasite communities and partially non-overlapping MHC diversity. While morphology explained 15% of the variation in mean parasite abundance, trophic ecology accounted for 16% and 22% of the MHC variation at the nucleotide and at the amino acid level, respectively. Parasitism and immunogenetic adaptation may thus add additional dimensions to the LT cichlid radiation.

Transgenerational plasticity and selection shape the adaptive potential of sticklebacks to salinity change.

In marine climate change research, salinity shifts have been widely overlooked. While widespread desalination effects are expected in higher latitudes, salinity is predicted to increase closer to the equator. We took advantage of the steep salinity gradient of the Baltic Sea as a space-for-time design to address effects of salinity change on populations. Additionally, genetic diversity, a prerequisite for adaptive responses, is reduced in Baltic compared to Atlantic populations. On the one hand, adaptive transgenerational plasticity (TGP) might buffer the effects of environmental change, which may be of particular importance under reduced genetic variation. On the other hand, physiological trade-offs due to environmental stress may hamper parental provisioning to offspring thereby intensifying the impact of climate change across generations (nonadaptive TGP). Here, we studied both hypothesis of adaptive and nonadaptive TGP in the three-spined stickleback (Gasterosteus aculeatus) fish model along the strong salinity gradient of the Baltic Sea in a space-for-time experiment. Each population tolerated desalination well, which was not altered by parental exposure to low salinity. Despite a common marine ancestor, populations locally adapted to low salinity lost their ability to cope with fully marine conditions, resulting in lower survival and reduced relative fitness. Negative transgenerational effects were evident in early life stages, but disappeared after selection via mortality occurred during the first 12-30 days posthatch. Modeling various strengths of selection, we showed that nonadaptive transgenerational plasticity accelerated evolution by increasing directional selection within the offspring generation. Qualitatively, when genetic diversity is large, we predict that such effects will facilitate rapid adaptation and population persistence, while below a certain threshold populations suffer a higher risk of local extinction. Overall, our results suggest that transgenerational plasticity and selection are not independent of each other and thereby highlight a current gap in TGP studies.

Disentangling Incomplete Lineage Sorting and Introgression to Refine Species-Tree Estimates for Lake Tanganyika Cichlid Fishes.

Adaptive radiation is thought to be responsible for the evolution of a great portion of the past and present diversity of life. Instances of adaptive radiation, characterized by the rapid emergence of an array of species as a consequence to their adaptation to distinct ecological niches, are important study systems in evolutionary biology. However, because of the rapid lineage formation in these groups, and occasional gene flow between the participating species, it is often difficult to reconstruct the phylogenetic history of species that underwent an adaptive radiation. In this study, we present a novel approach for species-tree estimation in rapidly diversifying lineages, where introgression is known to occur, and apply it to a multimarker data set containing up to 16 specimens per species for a set of 45 species of East African cichlid fishes (522 individuals in total), with a main focus on the cichlid species flock of Lake Tanganyika. We first identified, using age distributions of most recent common ancestors in individual gene trees, those lineages in our data set that show strong signatures of past introgression. This led us to formulate three hypotheses of introgression between different lineages of Tanganyika cichlids: the ancestor of Boulengerochromini (or of Boulengerochromini and Bathybatini) received genomic material from the derived H-lineage; the common ancestor of Cyprichromini and Perissodini experienced, in turn, introgression from Boulengerochromini and/or Bathybatini; and the Lake Tanganyika Haplochromini and closely related riverine lineages received genetic material from Cyphotilapiini. We then applied the multispecies coalescent model to estimate the species tree of Lake Tanganyika cichlids, but excluded the lineages involved in these introgression events, as the multispecies coalescent model does not incorporate introgression. This resulted in a robust species tree, in which the Lamprologini were placed as sister lineage to the H-lineage (including the Eretmodini), and we identify a series of rapid splitting events at the base of the H-lineage. Divergence ages estimated with the multispecies coalescent model were substantially younger than age estimates based on concatenation, and agree with the geological history of the Great Lakes of East Africa. Finally, we formally tested the three hypotheses of introgression using a likelihood framework, and find strong support for introgression between some of the cichlid tribes of Lake Tanganyika.

Back to Tanganyika: a case of recent trans-species-flock dispersal in East African haplochromine cichlid fishes.

The species flocks of cichlid fishes in the East African Great Lakes are the largest vertebrate adaptive radiations in the world and illustrious textbook examples of convergent evolution between independent species assemblages. Although recent studies suggest some degrees of genetic exchange between riverine taxa and the lake faunas, not a single cichlid species is known from Lakes Tanganyika, Malawi and Victoria that is derived from the radiation associated with another of these lakes. Here, we report the discovery of a haplochromine cichlid species in Lake Tanganyika, which belongs genetically to the species flock of haplochromines of the Lake Victoria region. The new species colonized Lake Tanganyika only recently, suggesting that faunal exchange across watersheds and, hence, between isolated ichthyofaunas, is more common than previously thought.

A tribal level phylogeny of Lake Tanganyika cichlid fishes based on a genomic multi-marker approach.

The species-flocks of cichlid fishes in the East African Great Lakes Victoria, Malawi and Tanganyika constitute the most diverse extant adaptive radiations in vertebrates. Lake Tanganyika, the oldest of the lakes, harbors the morphologically and genetically most diverse assemblage of cichlids and contains the highest number of endemic cichlid genera of all African lakes. Based on morphological grounds, the Tanganyikan cichlid species have been grouped into 12-16 distinct lineages, so-called tribes. While the monophyly of most of the tribes is well established, the phylogenetic relationships among the tribes remain largely elusive. Here, we present a new tribal level phylogenetic hypothesis for the cichlid fishes of Lake Tanganyika that is based on the so far largest set of nuclear markers and a total alignment length of close to 18kb. Using next-generation amplicon sequencing with the 454 pyrosequencing technology, we compiled a dataset consisting of 42 nuclear loci in 45 East African cichlid species, which we subjected to maximum likelihood and Bayesian inference phylogenetic analyses. We analyzed the entire concatenated dataset and each marker individually, and performed a Bayesian concordance analysis and gene tree discordance tests. Overall, we find strong support for a position of the Oreochromini, Boulengerochromini, Bathybatini and Trematocarini outside of a clade combining the substrate spawning Lamprologini and the mouthbrooding tribes of the 'H-lineage', which are both strongly supported to be monophyletic. The Eretmodini are firmly placed within the 'H-lineage', as sister-group to the most species-rich tribe of cichlids, the Haplochromini. The phylogenetic relationships at the base of the 'H-lineage' received less support, which is likely due to high speciation rates in the early phase of the radiation. Discordance among gene trees and marker sets further suggests the occurrence of past hybridization and/or incomplete lineage sorting in the cichlid fishes of Lake Tanganyika.

The evolution of cichlid fish egg-spots is linked with a cis-regulatory change.

The origin of novel phenotypic characters is a key component in organismal diversification; yet, the mechanisms underlying the emergence of such evolutionary novelties are largely unknown. Here we examine the origin of egg-spots, an evolutionary innovation of the most species-rich group of cichlids, the haplochromines, where these conspicuous male fin colour markings are involved in mating. Applying a combination of RNAseq, comparative genomics and functional experiments, we identify two novel pigmentation genes, fhl2a and fhl2b, and show that especially the more rapidly evolving b-paralog is associated with egg-spot formation. We further find that egg-spot bearing haplochromines, but not other cichlids, feature a transposable element in the cis-regulatory region of fhl2b. Using transgenic zebrafish, we finally demonstrate that this region shows specific enhancer activities in iridophores, a type of pigment cells found in egg-spots, suggesting that a cis-regulatory change is causally linked to the gain of expression in egg-spot bearing haplochromines.

A novel primer set for multilocus phylogenetic inference in East African cichlid fishes.

The cichlid fishes in the East African Great Lakes are a prime model system for the study of adaptive radiation. Therefore, the availability of an elaborate phylogenetic framework is an important prerequisite. Previous phylogenetic hypotheses on East African cichlids are mainly based on mitochondrial and/or fragment-based markers, and, to date, no taxon-rich phylogeny exists that is based on multilocus DNA sequence data. Here, we present the design of an extensive new primer set (24 nuclear makers) for East African cichlids that will be used for multilocus phylogenetic analyses in the future. The primers are designed to work for both Sanger sequencing and next-generation sequencing with the 454 technology. As a proof of principle, we validate these primers in a phylogenetically representative set of 16 cichlid species from Lake Tanganyika and main river systems in the area and provide a basic evaluation of the markers with respect to marker length and diversity indices.